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Abstract 

This  report  details  our  research  progress  during  the  course  of  this  grant.  This  report  includes 
development  of  new  wideband  antenna  elements,  strategies  for  cost-effective  random  arrays  using 
these  elements,  and  supporting  computational  codes.  First,  a  summary  of  the  development  and 
analysis  of  the  new  antenna  element  developed  under  this  project  is  provided.  Measured  and 
simulated  results  show  that  this  new  element  can  provide  a  3:1  bandwidth  (using  VSWR  =  3  as  the 
bandwidth  criteria  typically  used  for  wideband  elements)  with  no  degradation  in  radiation 
characteristics.  Next,  the  results  of  our  work  on  arrays  of  random  subarrays  are  presented. 
Theoretical  analyses  show  that  rotating  subarrays  composed  of  random  arrays  can  deliver 
performance  close  to  that  of  a  purely  random  array  while  making  the  array  more  cost  effective  and 
price  competitive  with  alternatives  that  work  over  much  smaller  bandwidths.  Parallel  work  on 
computational  tools  that  enable  wideband  simulation  of  elements  in  random  arrays  is  also  discussed. 
Conclusions  about  this  research  and  possible  directions  for  future  research  are  provided  along  with  a 
listing  of  the  publications  sponsored  by  this  project  to  date. 
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1.  Introduction 


Wideband  antennas  may  be  used  to  decrease  the  number  of  antennas  needed  on  platforms  that 
require  operation  over  a  large  frequency  band.  Some  typical  examples  of  wideband  antennas  include 
TEM  horns  [1,  2,  3]  and  Vivaldi  [4,  5,  6]  tapered  slot  antennas.  Although  both  types  do  offer 
wideband  operation,  they  are  often  not  suitable  for  certain  situations  due  to  their  large  thickness, 
weight,  and  complicated  feed  networks.  In  order  to  resolve  these  issues,  a  thin  low-profde  antenna 
is  needed.  The  attractive  features  of  canted  sector  antennas,  which  are  triangular,  conducting  sectors 
with  various  geometries  canted  above  ground  planes,  are  simplicity  and  wide  impedance  bandwidth. 
Band  limits  are  established  by  variations  in  the  radiation  pattern.  These  antennas  are  intended  for 
implementation  in  wideband  aperiodic  arrays  for  high-speed  communications. 

This  report  details  our  research  progress  during  the  course  of  this  grant.  This  work  includes 
development  of  new  wideband  antenna  elements,  strategies  for  cost-effective  random  arrays  using 
these  elements,  and  supporting  computational  codes.  Section  2  provides  a  summary  of  the 
development  and  analysis  of  the  new  antenna  element  developed  under  this  project.  Section  3 
summarizes  the  results  of  our  work  on  arrays  of  random  subarrays  that  promise  to  make  this 
architecture  cost  effective  and  price  competitive  with  less  functional  alternatives.  Parallel  work  on 
computational  tools  that  enable  wideband  simulation  of  elements  in  random  arrays  is  discussed  in 
Section  4.  Our  conclusions  about  this  research  and  possible  directions  for  future  research  are 
provided  in  Section  5  along  with  a  listing  of  the  publications  sponsored  by  this  project  to  date. 

2.  Wideband  Low  Profile  Antenna  Element 

The  project  began  with  examination  of  simple  canted  sector  antennas  [7,  8].  The  attractive  features 
of  canted  sector  antennas,  which  are  triangular,  conducting  sectors  with  various  geometries  canted 
above  ground  planes,  are  simplicity  and  wide  impedance  bandwidth.  Band  limits  are  established  by 
variations  in  the  radiation  pattern.  With  the  simple  sector,  the  radiation  pattern  exhibits  some 
depressions  close  to  broadside  that  are  unacceptable  for  adequate  scanning  of  the  array.  To  improve 
the  broadside  operation  of  the  element,  several  changes  were  implemented.  We  have  developed  a 
new  design  that  combines  a  number  of  geometrical  changes,  including  a  taper  in  the  main  triangular 
region  and  upward  bends  in  the  comers.  This  structure  is  shown  in  Figure  1  along  with  measured 
and  simulated  impedance  data.  From  2-6  GHz,  the  element  has  a  VSWR  that  is  always  less  than 
about  2.5:1. 
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Figure  1 :  Modified  canted  sector  antenna  (left)  with  measured  (blue-solid)  and  simulated  (orange-dotted) 

impedance  from  2-6  GFlz. 

Figure  2  shows  the  dimensions  of  the  tapered  sector  with  bent  corners.  The  piece-wise  taper  along 
the  side  edges  of  the  sector  antenna  is  defined  by  three  angles;  this  taper  resembles  a  quadratic 
function.  The  corners  were  bent  upwards  by  an  angle  C,  with  respect  to  the  plane  of  the  sector,  and  a 
cant  angle  of  a  =  5°  was  used.  The  same  ground  plane  and  excitation  method  used  for  the  simple 
sector  was  applied. 


Figure.  2:  Geometry  of  the  tapered  sector  with  bent  comers. 

The  antenna  was  first  simulated  with  corner  bends  of  C,  =  0°  (no  bend),  20°,  and  40°.  As  expected,  the 
input  impedance  improved  with  increasing  The  case  with  C,  =  20°  had  the  best  match  of  the 
simulated  cases.  The  simulated  radiation  patterns  showed  the  removal  of  the  depressions  at  3.6  and 
4.4  GHz  with  >  20°.  However,  the  comer  bends  were  not  able  to  remove  the  null  at  2  GHz.  Even 
so,  there  is  good  broadside  behavior  from  2.2  GHz  to  at  least  6  GHz;  simulations  above  6  GHz  were 
not  performed.  As  expected,  the  antenna  continues  to  exhibit  good  impedance  characteristics  as  the 
frequency  increases.  A  corner  bend  angle  of  C,  =  20°  was  sufficient  in  the  simulations  to  satisfy  the 
broadside  pattern  conditions,  and  maintain  the  low-profile  property  (this  antenna  has  a  vertical 
dimension  of  22.3  mm  and  the  simple  sector  with  a  =  10°,  /?  =  67°  and  L  =  125  mm  in  [9]  has  a 
vertical  dimension  of  22.1  mm),  and  was  thus  fabricated  and  is  shown  in  Figure  3.  The  element 
taper  reduces  the  VSWR  at  the  lower  end  of  the  band  compared  to  that  of  the  simple  sector. 
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Figure  3:  Fabricated  tapered  sector  with  bent  comers. 


Figure  4:  VSWR  of  the  tapered  sector  with  bent  comers. 


E-plane  H-plane 

0°  0° 


Figure  5:  Measured  and  simulated  (FGP)  radiation  patterns  at  3.2  GHz  of  the  tapered  sector  with  bent  comers. 
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E-plane  H-plane 

0°  0° 


Figure  6:  Measured  and  simulated  (FGP)  radiation  patterns  at  4.4  GFlz  of  the  tapered  sector  with  bent  comers. 

Figure  4  presents  the  measured  and  simulated  VSWR.  There  was  a  slight  over-estimation  in  the 
simulated  results  compared  to  the  measured  results,  but  the  overall  input  impedance  trends  agree 
well  with  each  other.  Figures  5  and  6  show  the  measured  and  simulated  (E-plane  copolar  only) 
patterns  at  3.2  GHz  and  4.4  GHz,  respectively,  with  the  FGP  designation  indicating  simulated  results 
obtained  with  the  TD-AIM  simulator  developed  in  this  project.  The  depressions  that  were  present  in 
the  simple  sector  were  indeed  alleviated  and  a  continuous  broadside  pattern  bandwidth  with  a  2.5:1 
VSWR  from  2.2-6  GHz  was  achieved. 

The  antenna  has  been  analyzed  extensively.  We  are  presently  finishing  an  analytical  model  of  both 
the  simple  sector  as  well  as  this  most  promising  permutation  to  support  new  designs  and  scaling  in 
the  future.  Two  journal  papers  on  this  element  are  in  preparation.  In  the  course  of  this  research,  Mr. 
Garvin  Cung,  the  graduate  student  who  developed  this  antenna  element,  completed  his  M.S.  thesis 
on  this  topic  [10]  and  graduated  from  the  University  of  Illinois  at  Urbana-Champaign  in  August 
2005.  He  is  presently  employed  with  Northrop  Grumman  in  Baltimore,  Maryland. 


3.  Random  Array  Studies 

While  truly  random  arrays  may  prove  too  costly  to  fabricate  efficiently  using  today’s  methods,  we 
are  developing  design  approaches  that  provide  good  sidelobe  performance  over  frequency  while 
remaining  practical  and  cost-effective  to  fabricate  and  operate.  Our  approach  is  to  use  a  smaller 
random  array  as  a  subarray  building  block.  A  brief  discussion  of  the  results  of  three  different 
variations  is  provided  here. 

We  have  studied  periodic  arrays  of  random  subarrays  (PARS)  that  consist  of  small  rectangular 
random  subarrays  arranged  periodically  [11].  Such  a  configuration  is  shown  in  Figure  7(a).  For  an 
array  of  a  given  aperture  size  and  number  of  elements,  the  robust  wideband  performance  of  a 
random  array  is  preserved  in  a  geometrically  simpler  design,  but  with  somewhat  higher  sidelobes 
than  a  purely  random  array  of  the  same  size.  As  a  result,  from  a  performance  standpoint  a  purely 
random  array  is  more  effective,  and  the  reason  to  use  a  PARS  would  be  that  one  wanted  to  use  a 
certain  number  of  elements  to  achieve  some  desired  gain  in  a  geometrically  simpler  package  than  the 
purely  random  array. 
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Another  approach  is  to  rotate  the  subarrays  relative  to  one  another  in  the  array,  creating  arrays  of 
periodically  rotated  random  subarrays  (ARRS-P),  depicted  in  Figure  7(b)  [12,  13].  Our  work  shows 
that  if  the  number  of  elements  in  the  subarrays  is  larger  than  a  threshold  value  for  a  particular 
frequency  range,  both  the  PARS  and  ARRS  have  very  stable  array  factors  and  can  easily  be  designed 
so  that  their  sidelobe  levels  remain  below  a  specified  value  over  the  entire  operating  band. 
However,  the  ARRS  has  a  lower  threshold  value,  as  shown  in  Figure  8,  resulting  in  lower  element 
density  and  consequently  an  extended  lower  frequency  limit.  The  number  of  elements  required  to 
achieve  a  specified  sidelobe  level  is  reduced,  which  results  in  arrays  with  fewer  elements  and  feed 
lines  for  the  same  sidelobe  level  or  in  lower  sidelobe  levels  for  the  same  number  of  elements.  Using 
fewer  elements  will  result  in  lower  aperture  gain,  but  may  also  result  in  a  cost  savings  over  multiple 
traditional  periodic  arrays  that  perform  comparable  functions. 
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Figure  7:  (a)  Periodic  array  of  random  subarray  (PARS),  (b)  Array  of  rotated  random  subarrays  (  ARRS). 
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Figure  8:  Number  of  elements  required  to  achieve  a  -20  dB  sidelobe  level  with  85%  probability  vs.  a 
(subarray  side  length  in  wavelengths),  which  depends  linearly  on  frequency,  for  PARS  (dashed  line)  and  an 

ARRS-P  (solid  line)  with  four  subarrays  each  [12]. 

Finally,  the  concept  of  subarray  rotation  is  generalized  to  allow  random  rotations.  In  an  array  of 
randomly  rotated  random  subarrays  (ARRS-R),  the  pattern  of  rotations  is  not  deterministic.  For 
simplicity,  we  require  that  the  subarrays  not  overlap.  As  with  the  PARS  and  ARRS-P,  the  likely 
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behavior  of  the  arrays  can  be  predicted  using  probabilistic  methods.  Results  indicate  that  the  ARRS- 
R  can  achieve  lower  sidelobe  levels  with  fewer  elements  per  subarray  as  more  subarrays  are  added 
[14,  15].  This  is  in  contrast  to  the  PARS  and  the  ARRS-P,  in  which  addition  of  more  subarrays 
eventually  ceases  to  improve  the  sidelobe  level  perfonnance  of  the  complete  array. 

A  journal  paper  on  this  random  array  study  has  been  published  [14]  and  is  attached  to  this  report  as 
Appendix  1.  In  the  course  of  this  research,  Ms.  Kiersten  Kerby  completed  her  M.S.  thesis  on  this 
topic  [15]  and  earned  the  M.S.  degree  from  the  University  of  Illinois  at  Urbana-Champaign  in  May 
2005.  She  spent  the  summer  of  2005  in  an  engineering  internship  at  Ball  Aerospace  in  Colorado, 
and  she  is  now  continuing  for  her  Ph.D.  at  the  University  of  Illinois  under  the  direction  of  Prof. 
Bernhard. 


We  requested  a  no-cost  extension  in  order  to  complete  a  measurement  demonstration  of  these 
random  subarray  concepts  using  optical  techniques.  The  screen  designs  were  generated  using  our 
established  procedures,  and  slots  are  used  as  the  radiating  elements.  Due  to  unforeseeable  delays 
related  to  the  breakdown  of  fabrication  equipment,  we  were  delayed  in  fabricating  the  first  designs 
until  May  2006.  After  fabrication,  we  determined  that  the  conductive  coating  on  the  screens  was  not 
thick  enough,  allowing  light  to  pass  through.  We  are  currently  working  on  a  slightly  modified 
fabrication  technique  that  will  result  in  thicker  conducting  film  deposition  in  the  optical  screens. 
Once  the  new  screens  are  fabricated,  we  will  measure  the  beam  patterns  produced  using  uniform 
laser  illumination  to  assess  the  feasibility  of  several  of  our  proposed  approaches  for  arrangement  of 
random  subarrays  in  large  arrays. 

4.  Computational  Tool  Development 

A  time  domain  integral  equation  solver  that  is  being  used  to  quickly  and  accurately  simulate  these 
antennas  has  been  developed  at  the  University  of  Illinois  under  this  grant.  Specifically,  a  marching- 
on-time  (MOT)  algorithm  for  solving  a  time-domain  electric  field  integral  equation  pertinent  to  the 
analysis  of  free-standing  antennas  was  constructed.  Given  a  bandlimited  excitation,  the  algorithm 
solves  for  induced  surface  currents  on  the  conductors  using  a  time  domain  method-of-moments 
formulation.  For  a  geometry  modeled  using  Ns  surface  unknowns,  the  computational  complexity  of 

this  algorithm  scales  as  0(NtN^)  for  Nt  time  steps  of  simulation.  Using  a  parallel  FFT-based 
accelerator,  the  time-domain  adaptive  integral  method  (TD-AIM),  the  cost  of  solution  is  reduced  to 
0(NtNs  log2  Ns)  operations  for  the  antennas  shown  below.  The  solver  employs  the  message  passing 
paradigm  via  the  MPI  standard  for  communication  between  processors  and  utilizes  the  Fastest 
Fourier  transform  in  the  West  (FFTW)  library  to  compute  parallel  FFTs.  Thanks  to  the  near-ideal 
parallel  performance  of  the  solver,  we  have  been  able  to  analyze  the  broadband  characteristics  of 
antennas  that  involve  Ns  >100,000  degrees  of  freedom  using  supercomputers  with  tens  of  processors 
in  practical  amounts  of  time. 

Various  antenna  shapes  have  been  investigated  using  the  TD-AIM  based  simulator  to  investigate 
their  radiation  patterns  over  a  broad  range  of  frequencies  [16,  17].  The  simulations  shown  in  Section 
2  with  finite  ground  planes  were  performed  with  this  simulator  and  agree  well  with  measured  results. 


Further  array  simulations  will  rely  mainly  on  array  factor  multiplications  of  element  patterns  due  to 
the  otherwise  extremely  large  computational  space  required. 

5.  Conclusions  and  Directions  for  Future  Work 

Our  work  on  the  canted  sector  antenna  led  to  modified  designs  with  good  impedance  and  radiation 
characteristics  that  can  meet  a  number  of  different  array  scan  requirements.  We  continue  to  work  on 
analytical  models  for  all  of  these  designs  to  enable  further  refinement  as  well  as  scaling  for  new 
frequency  bands  of  interest.  In  the  future,  more  rigorous  studies  on  the  corner  bend  locations  and 
different  quadratic  tapers  will  be  pursued  to  improve  the  antenna  perfonnance  further. 

Our  work  in  random  arrays  has  shown  very  promising  results.  Work  on  the  possible  packaging  and 
arrangement  of  random  subarrays  may  continue,  including  investigation  of  the  effects  of  small 
numbers  of  elements  and  the  possibility  of  recursive  rotation  arrangements  that  may  duplicate  the 
behavior  of  truly  random  arrays. 

Finally,  we  will  continue  to  pursue  the  measurement  demonstration  using  optical  techniques  during 
the  last  period  of  this  project  we  will  calculate  complete  radiation  patterns  for  arrays  with  random 
subarrays  composed  of  these  wideband  elements.  During  this  process,  we  will  also  evaluate  scan 
capabilities  and  investigate  the  possibility  of  using  the  wideband  element’s  endfire  characteristics  at 
higher  frequencies  to  extend  scan  angle  ranges. 
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Abstract — We  propose  an  extension  to  the  concept  of  random 
antenna  arrays.  Three  new  array  geometries  are  introduced 
which  are  intended  to  incorporate  the  wide  bandwidth  capability 
of  a  random  array  while  having  simplified  geometry  to  increase 
their  suitability  for  lower-cost  applications.  The  behaviors  of  the 
periodic  array  of  random  subarrays,  array  of  periodically  rotated 
random  subarrays,  and  array  of  randomly  rotated  random  subar¬ 
rays  are  each  characterized  probabilistically.  Results  indicate  that 
subarray  rotation  can  lower  the  sidelobe  level  of  the  array  factor. 
An  example  of  the  design  process  is  presented  with  calculated 
array  factors. 

Index  Terms — Antenna  arrays,  random  arrays,  subarrays,  wide¬ 
band  arrays. 


I.  Introduction 

FOR  many  applications,  it  is  common  to  use  periodic  an¬ 
tenna  arrays  because  they  are  relatively  simple  to  design 
and  manufacture.  However,  as  the  importance  of  wide  band- 
widths  and  multifrequency  operation  has  increased,  it  has  be¬ 
come  evident  that  there  are  significant  disadvantages  to  the  use 
of  periodic  arrays.  For  instance,  because  the  element  spacing  of 
periodic  arrays  is  proportional  to  the  wavelength  of  operation, 
the  patterns  of  periodic  arrays  change  drastically  and  quickly  de¬ 
velop  grating  lobes  in  response  to  changes  in  frequency.  In  the 
past,  these  have  not  been  major  issues  because  systems  typically 
operated  at  only  one  frequency  and  with  narrower  band  widths. 
But  as  the  trend  of  increasing  bandwidths  and  multifrequency 
systems  continues,  it  will  be  important  for  array  design  to  over¬ 
come  those  obstacles.  Aperiodic  element  placement  offers  one 
possible  solution. 

When  aperiodic  arrays  were  first  investigated  in  the  1960s, 
their  utility  was  based  largely  in  their  capability  for  high  reso¬ 
lution  with  relatively  few  elements,  as  well  as  wide  scan  angle 
capability.  They  were  considered  a  lower-cost  replacement 
for  the  large,  mechanically  scanned  reflector  antennas  used  in 
remote  sensing.  Initially,  designers  of  aperiodic  arrays  often 
used  preassigned  element  locations,  in  essentially  arbitrary 
distributions  [1],  Because  of  the  limitations  of  computers  at  the 
time,  optimization  was  somewhat  impractical  and  attempts  at  it 
led  to  only  locally  optimum  element  placement  [2],  [3],  Lo  [4] 
then  proposed  treating  large  antenna  arrays  probabilistically, 
even  though  the  problem  was  not  originally  probabilistic.  This 
approach — probabilistically  analyzing  arrays  with  random 
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element  placement — represented  the  formal  introduction  of 
random  arrays  and  a  systematic  method  of  analyzing  them.  In¬ 
vestigation  of  random  arrays  produced  a  general  understanding 
of  their  behavior  and  advantages  [5]  and  found  that,  in  fact, 
most  deterministic  aperiodic  array  design  methods  produced 
results  that  were  not  significantly  better  than  random  [6]. 
Arrays  that  were  quasi-optimized  using  dynamic  programming 
were  the  only  ones  whose  performance  was  appreciably  better 
than  random  element  placements.  Further  work  on  the  subject 
included  some  interesting  special  cases  of  the  purely  random 
array,  such  as  binned  [7]  and  nearest-neighbor  constrained  [8] 
random  arrays. 

Aperiodic  arrays  have  very  stable  patterns  over  wide  ranges 
of  frequency,  so  they  are  well  suited  to  operating  over  large 
bandwidths.  However,  the  aperiodicity  that  yields  this  behavior 
also  makes  such  arrays  difficult  to  fabricate.  It  is  far  simpler 
to  manufacture  a  repetitive  structure  such  as  a  periodic  array. 
Recently,  work  on  aperiodic  arrays  has  often  been  concerned 
with  optimization  of  their  patterns.  Naturalistic  optimization 
methods  such  as  genetic  [9]  and  particle-swarm  [10]  algorithms 
are  certainly  capable  of  producing  arrays  with  performance 
better  tailored  to  their  particular  applications.  However,  since 
the  problem  considered  here  is  largely  one  of  subarray  arrange¬ 
ment,  this  work  addresses  random  arrays.  Once  some  effective 
geometric  simplification  methods  are  found,  they  may  later  be 
combined  with  optimization  techniques.  Since  random  arrays 
by  definition  choose  one  element  distribution  out  of  the  entire 
set  of  all  possible  arrays,  this  approach  can  be  applied  to  other 
types  of  aperiodic  arrays. 

The  following  sections  detail  attempts  at  geometry  simplifi¬ 
cation,  including  a  periodic  array  of  random  subarrays  (PARS) 
and  arrays  of  periodically  and  randomly  rotated  random  subar¬ 
rays  (ARRS-P  and  ARRS-R,  respectively).  A  comparison  be¬ 
tween  the  three  types  is  presented,  followed  by  remarks  on  the 
results  of  the  investigation  and  possibilities  for  future  work. 

II.  Three  Geometries  for  an  Array  of  Random 
Subarrays 

The  random  element  placement  of  purely  random  arrays 
makes  them  more  difficult  to  manufacture  than  a  periodic 
array.  Three  array  geometries  were  developed  in  an  attempt 
to  mitigate  this  problem  while  preserving  the  wide  bandwidth 
response  and  high  failure  tolerance  of  purely  random  arrays. 
This  section  presents  the  geometries  of  all  three  types. 

The  first  type  of  array,  a  periodic  array  of  random  subarrays 
(PARS),  is  composed  of  many  identical  random  subarrays  ar¬ 
ranged  in  a  periodic  grid  (Fig.  1). 
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Fig.  1.  Geometry  of  PARS. 
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Fig.  2.  ARRS-P  geometry. 

The  next  variation  is  an  array  of  periodically  rotated  random 
subarrays  (ARRS-P).  The  identical  subarrays  of  the  PARS 
are  rotated  in  the  azimuthal  plane  in  a  pattern  that  repeats 
over  blocks  of  four  subarrays.  In  each  block,  the  coordi¬ 
nates  of  the  subarrays’  centers  are  (—(a/2),  a/2),  (a/2,  a/2), 
((a/2), -(a/2)),  and  (—(a/2), -(a/2)),  respectively.  The 
subarrays  are  rotated  in  </  by  0,  90,  180,  and  270°,  respec¬ 
tively.1  It  is  assumed  that  the  array  elements  have  rotationally 
symmetric  radiation  patterns.  Fig.  2  shows  one  such  block. 

Last,  the  concept  of  subarray  rotation  is  generalized  to  allow 
random  rotations.  In  an  array  of  randomly  rotated  random  sub¬ 
arrays  (ARRS-R),  the  pattern  of  rotations  is  not  deterministic. 
The  allowed  rotation  values  may  be  discrete — for  example,  in 
Fig.  3  the  rotations  allowed  are  0,  90,  180,  and  270° — or  they 
may  be  a  continuum.  For  this  analysis  we  will  require  that  the 
subarrays  not  overlap. 

For  all  three  arrays,  a  and  b  are  subarray  dimensions,  mea¬ 
sured  in  wavelengths,  and  P  and  Q  are  the  numbers  of  subar¬ 
rays  in  the  <j>  =  0°  and  <u  =  90°  directions,  respectively.  For 
the  rectangular  subarrays  of  the  PARS,  a  and  b  denote  the  sides 
of  the  subarrays.  Due  to  its  rotations,  the  ARRS-P  is  required  to 

1  When  individual  circularly  polarized  elements  are  rotated  in  the  pattern 
shown  for  the  ARRS-P,  grating  lobes  arise  due  to  the  introduced  periodicity. 
The  same  rotation  pattern  was  applied  again  to  the  four-element  subarrays,  and 
it  was  found  that  when  four  subarrays  were  arranged  in  a  square  sixteen-ele¬ 
ment  array,  the  grating  lobes  were  reduced  [11].  This  work  employs  the  same 
rotation  pattern  in  order  to  demonstrate  the  effects  of  non-periodic  (a  single 
group  of  four  subarrays)  and  periodic  (more  than  one  group)  rotation. 
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Fig.  3.  Geometry  of  a  particular  ARRS-R  array,  with  round  subarrays  and 
M  =  4. 


Fig.  4.  Illustration  of  variable  definitions.  Left:  subarray  indices  and  global 
coordinate  system.  Right:  subarray  rotation,  dimensions,  and  local  coordinate 
system. 


have  square  subarrays  whose  side  length  is  denoted  a.  In  con¬ 
trast,  the  ARRS-R  subarrays  may  take  any  shape,  so  a  is  de¬ 
fined  to  be  the  subarrays’  largest  dimension.  N  is  the  number 
of  elements  per  subarray.  The  random  variables  xn  and  yn  de¬ 
note  the  normalized  positions  of  individual  elements  in  the  sub¬ 
array  in  a  coordinate  system  local  to  the  subarray.  They  may 
vary  from  —1/2  to  1/2.  xn  and  yn  have  probability  distributions 
f(x )  and  g(y),  respectively.  Both  distributions  shall  be  required 
to  be  even,  because  it  reduces  the  complexity  of  calculations  and 
because  most  arrays  in  practice  have  symmetric  element  place¬ 
ments.  In  the  ARRS-P  and  ARRS-R,  the  angle  by  which  a  sub¬ 
array  with  indices  p,  q  is  rotated  is  denoted  ippq.  The  definitions 
of  these  variables  are  illustrated  in  Fig.  4.  For  the  ARRS-R,  we 
must  also  define  M,  the  number  of  allowed  values  for  the  rota¬ 
tion  angle. 

III.  Array  Factor  Analysis 
A.  Array  Factors 

For  each  array  geometry,  the  array  factor  is  the  superposition 
of  the  patterns  of  individual  isotropic  elements.  For  the  PARS, 
this  immediately  simplifies  into  the  familiar  sine  function  of  the 
periodic  superarray  multiplied  by  the  sum  associated  with  the 
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subarrays,  as  is  apparent  in  this  expression  for  the  array  factor 
GpARS  <f>) 


part  is  then  written 


Gpars  = 


1  sin(P7rasin#cos</>) 
NPQ  sin(7ra  sin  6  cos  <j>) 


G  ARRS-P, #($,0) 


x  All(Q7rb  sill  6  sill  (l>)  \  gj27T  sin  0(axn  cos  4)Jrbyri  sin  d)  f  |  j 

sin(-7r&sin0sinfl!))  ^ 

v  T 7  n=l 

To  begin  with,  we  will  consider  a  single  four-subarray  block 
of  the  ARRS-P,  and  later  generalize  to  arrays  of  many  subarrays. 
For  a  single  block,  the  array  factor  is  given  by  the  sum  of  the 
subarray  factors,  with  phases  adjusted  for  rotation  and  position. 
The  array  factors  of  subarrays  rotated  by  0,  90,  180,  and  270° 
are  denoted  Go,  Ggo,  Giso,  and  G270,  respectively 


_  ___gj7rct  sin  0( —  cos  0-)-sin  <$>) 

X  E  gj  2tt a  sin  6(xn  cos  4>-\-yn  sin  <f>) 
n 

q  _  jj_  gj7ra  sin  6(c os  0+sin  <^>) 


cos  (-7T a  sin  0(—  cos  <j>  +  sin  <f>  +  2.?;r,  cos  (p 

n= 1 

+2t/ra  sin  0)) 

N 

+7^  cos  (7ra  sin  #(cos  <f>  +  sin  <j>  —  2xn  sin  <>> 


+^yn  cos  <f>))  J .  (8) 

Finally,  the  real  and  imaginary  parts  of  the  ARRS-R  array 
factor  are  given  in  (9)  and  (10)  as  follows: 


Garrs-r,r 


p  Q  N 

EEE  cos(27ra  sin  8(p  cos  <f>  +  q  sin  < 


NPQ  ^  ^  ^ 

p= 1 q= 1 n=l 


x  ej27ra  sin  sin  cos  (3)  Garrs-r,/ 


H“  C0S((^  'Ippq) 

+  y„  sin -  V’pg))  (9) 


P  Q  N 


G180  =  —ej^asin0(.cos<t—sin<j>) 


x  ef27ra  sin  e(-a;„  cos  <j>—yn  sin  <j>)  ^ 


Gty-Q  ^  ^,77m  sin  —  cos  d — sin  <: 


X  E  gj2na  sin  6(xn  sin  4>—yn  cos  <^>) 


GaRRS-p(^,  4)  —  ^(Go  +  Ggo  +  Giso  +  G27o)-  (6) 

The  array  factor  of  the  ARRS-R  is  similarly  constructed  from 
a  sum  of  subarray  factors 


G( A  a\  _  ^  \  ^  \  A  ^7'27rasin0(pcos<£+<7sin</>) 

ARRS— R(t 'iV)  ~  JypQ  2^l^e 
*  P= 1  9=1 
N 

x  gj‘27ra  sin  ««(o— yr,,)  •  siti(.;— i/:,,,,))  ^ 


=  w?oEEE  sin(27ra  sin  #(p  cos  <j>  +  g  sin  4> 

^  p=l  q=l  n=l 

+  :i:n  cos(0  -  |!ipq) 

T  Vn  sin(</>  t/ipg)).  (10) 

To  compute  the  expected  value  and  variance  of  these  quanti¬ 
ties,  we  will  introduce  reduced  angle  variables  u  and  t>,  and  char¬ 
acteristic  functions  ipx(u)  and  <py(u).  The  characteristic  func¬ 
tions  are  only  coincidentally  similar  to  the  functions  used  in  sta¬ 
tistical  analysis;  they  actually  arise  because  they  are  related  to 
the  expected  array  factor  of  a  subarray 

u  =  sin  6  cos  4>  (11) 

v  =  sm6sm(j)  (12) 

1/2  1 

ip:E(u)  =  J  f(x)e?2’*auxdx  =  j  f(x)cos(2iraux)dx  (13) 

-1/2  -l 

1/2  i 

<Py(v)=  J  g(y)ej2”bvy dy  =  j  g(y)cas(2irbvy)dy.  (14) 


B.  Expected  Value  and  Variance 

Although  the  element  positions  are  as  yet  unknown,  the 
likely  behavior  of  the  arrays  can  be  predicted  using  proba¬ 
bilistic  methods  [4],  In  this  analysis,  we  will  characterize  the 
approximate  shape  of  the  array  factor.  It  is  most  useful  to  treat 
the  real  and  imaginary  parts  of  each  array  factor  separately.  The 
real  part  of  the  PARS  array  factor  will  be  denoted  Gpars, R, 
the  imaginary  part  Gpars, J,  and  similar  symbols  are  assigned 
to  the  real  and  imaginary  parts  of  the  other  arrays’  factors. 

In  the  case  of  the  ARRS-P,  the  particular  rotation  pattern 
chosen  causes  the  imaginary  part  to  be  identically  zero.  The  real 


Since  f(x)  and  g(y)  are  even  functions,  the  integral  of  the  ex¬ 
ponential  is  equivalent  to  the  integral  of  the  cosine. 

Because  the  ARRS-R  is  more  general,  it  will  require  slightly 
modified  expressions.  Since  the  cosine  and  sine  factors  on  x  and 
y  occur  in  pairs,  we  will  define  a  reduced  angle  variable  sum  Um 
to  replace  u  and  v.  Below,  ipm  is  the  mth  allowed  rotation  angle 
out  of  the  total  M 

Um  =  sin (9  (xcos(<f)  -  ?/>m)  +  ysin(</>  -  $TO))  (15) 
=  JJcos(2,aUm)f(x)g(y)dA.  (16) 


2256 


IEEE  TRANSACTIONS  ON  ANTENNAS  AND  PROPAGATION,  VOL.  54,  NO.  8,  AUGUST  2006 


The  quantity  (p(Um(x,  y))  above  is  similar  to  the  ip  definitions 
used  earlier,  but  the  expression  here  is  more  general  and  does  not 
separate  the  x  and  y  integrals,  since  the  cosine  and  sine  factors 
on  x  and  y  occur  in  pairs  that  depend  on  the  value  of  p  for  the 
subarray. 

With  these  definitions  in  mind,  the  array  factor  expected 
values  are  as  follows.  For  the  PARS 


E  (Gpars,r(u,  v)) 


1  sm(iraPu )  sin(67r(3u) 
PQ  sin(7rart)  sin(&7ru) 
X  <px(2TTau)(py(2nav) 


E  (Gpars,i(u,v))  =0. 


(17) 

(18) 


For  the  ARRS-P 


E  (G arrs-p, r(u,  v))  =  -  (cos  (?r a(u  -  v))  ipx(2trau) 
xipy(27rav)  +  cos(na(u  +  v))  ipx(27rav)ipy(27rau)) .  (19) 

And  for  the  ARRS-R 


1 

PQM 


E(Garrs-r,r) 

( P,Q )  M 

E  E  cos  (27ra  (pu+qv+Um(  x,  y))) 

A  (I’-'/MM)"'  1 

X  f(x)g(y)dA  (20) 

E(Garrs-r,i) 

1  ,  ,  (PQ)  m 

=  7 Am  E  E  sin  (2-7T a  ( pu+qv+Um(x ,  y))) 

4  J  A  (P-9)=( hi)rn=l 

x  f(x)g(y)dA.  (21) 


In  (20)  and  (21),  A  is  the  aperture  of  a  subarray  and  dA  is  a 
differential  area  element.  If  the  shape  of  the  subarrays  is  sym¬ 
metric  in  x  and  y,  then  dA  is  an  even  function  of  both,  and  these 
expressions  simplify  to 


1  (P,Q)  M 

E(G arrs-r, r)  =  pn ,  r  E  E  cos  (2ira(pu  +  qv)) 

V  (p,g)=(l,l)  m=l 

x  c p(2'KaUm )  (22) 

^  (P,Q)  M 

E(Garrs-rj)  =  pqm  E  E  sin  (2na(pu  +  qv)) 

W  (jp,q)=(  1,1)  ™=1 

x  <p(2'KaUm).  (23) 

In  the  region  of  the  main  beam,  the  array  factor  is  nearly  de¬ 
terministic.  It  is  very  likely  to  have  the  shape  of  the  array  factor’ s 
expected  value  [12].  In  the  sidelobe  region  this  is  not  the  case, 
so  it  remains  for  us  to  characterize  the  sidelobe  level  (SLL)  in 
order  to  arrive  at  an  accurate  idea  of  the  array’s  behavior.  It  is 
helpful  that,  by  the  Riemann-Lebesque  theorem,  <p  terms  will 
go  to  zero  in  the  sidelobe  region.  As  a  result,  in  that  region  all 
of  the  expected  values  are  approximately  zero. 

By  the  same  reasoning,  approximations  can  also  be  made  to 
the  variances  of  the  real  and  imaginary  parts  of  each  array  factor 


in  the  sidelobe  region.  The  approximate  variances  are  given 
below.  For  the  PARS 

2  _  2 
^PARS,i?  — ^PARSJ 

/  1  sin(a7rPu)  sin(67rQiQ \2  1 

\PQ  sin(a7 tu)  sm(bnv)  )  2 N‘ 

Since  the  imaginary  part  of  the  ARRS-P  array  factor  is  identi¬ 
cally  zero,  the  variance  of  the  imaginary  part  must  be  zero.  The 
variance  of  the  real  part  is 


a 


1 

AN' 


(25) 


Lastly,  the  variances  for  the  real  and  imaginary  parts  of  the 
ARRS-R  array  factor  are 


2  _  2 
°ARRS— R,i?  —  °ARRS-R,J 

1  1 

~  2 NPQ  +  2 NMP2Q2 

sin  ( airu(P  +  1)) 

sin(a7ra) 

sin  ( aTTv(Q  +  1))  \  2 
sin(a7ru)  J 

In  the  special  case  where  M  is  very  large  (corresponding  to 
a  continuum  of  allowed  rotation  angles),  the  variances  for  the 
ARRS-R  are  approximately  equal  and  given  by  the  simpler  ex¬ 
pression 


PQ  ■  (26) 


ctARRS-R,R  ~  ^ARRS-R,/ 


1 

2  NPQ' 


(27) 


C.  Probability  Distributions 

By  the  central  limit  theorem,  if  N  (and,  in  the  case  of  the 
ARRS-R,  the  product  PQ)  is  large,  then  each  of  the  real  and 
imaginary  parts  of  all  of  these  array  factors  is  an  approximately 
normally  distributed  random  quantity.  In  that  case,  the  magni¬ 
tudes  of  the  PARS  and  ARRS-R  array  factors  at  any  particular 
observation  angle  have  chi-squared  distributions  with  two  de¬ 
grees  of  freedom,  and  the  sidelobes  of  the  ARRS-P  array  factor 
have  a  normal  distribution.  This  can  be  used  to  arrive  at  a  prob¬ 
ability  distribution  for  the  height  of  a  particular  sidelobe. 

For  the  PARS,  the  probability  that  a  sidelobe  is  less  than  some 
fraction  r  of  the  main  beam  is  given  by 


Pr{|GpARs(tii,'6’i)|  <r}  =  1  -  e  r'/^\  (28) 

Since  ctp  =  cry,  we  can  use  a  single  quantity  a2. 

In  the  case  of  the  ARRS-P,  the  probability  distribution  of  the 
sidelobe  magnitude  is 


Pr{|GARRs-p(tii5:ti)|  <  r}  ss  Erf (rV2N).  (29) 
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When  considering  the  ARRS-R,  we  may  either  look  at  the 
case  where  M  is  large  and  we  can  use  the  approximate  vari¬ 
ance  from  (27),  or  the  case  where  M  is  not  large  and  we  must 
use  the  variance  from  (26).  In  terms  of  utility,  we  are  more  in¬ 
terested  in  the  case  where  M  is  large  because  in  that  case  the 
whole  variance  depends  on  1/PQ.  If  M  is  small  only  one  of  the 
terms  depends  on  1/PQ — the  sine 2  in  the  second  term  some¬ 
times  cancels  out  the  P2Q 2  in  the  denominator,  depending  on 
the  angle.  Because  of  this,  the  case  of  large  M  is  the  more  inter¬ 
esting  case  and  from  now  on  it  will  be  the  only  type  of  ARRS-R 
considered.  The  probability  distribution  of  the  magnitude  of  a 
sidelobe  for  this  case  is 


Pt{|G(«i,vl)|  <  r}  =  1  -  e-NPQr\  (30) 

From  these  probability  distributions,  the  expected  value  of 
the  SLL,  which  predicts  the  average  SLL,  can  be  computed  for 
each  type  of  array.  Since  the  PARS  SLL  has  a  chi-squared  dis¬ 
tribution,  it  is  known  that  the  average  SLL  for  a  given  location 
must  be  2<tpars-  However,  since  the  variance  of  the  PARS  SLL 
is  angle-dependent,  we  must  also  average  over  the  solid  angle 
occupied  by  the  sidelobes.  Assuming  the  main  beam  is  narrow, 
we  can  express  the  PARS  average  SLL  as  follows: 

ravg  =  ^  J  J (2(Jpars)  sin  OdOdcj)  (3 1 ) 

The  sidelobe  value  of  the  ARRS-P  has  a  normal  distribution,  and 
we  must  take  its  absolute  value  in  order  to  calculate  the  average 
SLL  in  a  useful  way 


—  OO 


The  average  SLL  of  the  ARRS-R  is  obtained  directly  from  the 
properties  of  a  chi-squared  distribution  with  two  degrees  of 
freedom: 


The  probability  distribution  of  the  peak  SLL  can  be  approxi¬ 
mately  characterized  by  a  joint  probability  distribution  of  many 
individual  sidelobes  [12].  Because  of  the  angle-dependent  vari¬ 
ance  of  the  PARS,  it  is  clear  that  the  highest  sidelobes  are  likely 
to  be  at  the  locations  of  grating  lobes  in  the  superarray  pattern. 
At  those  locations,  the  variance  is  1  /2 N  and  the  probability  dis¬ 
tribution  of  the  PARS  peak  sidelobe  level  is 


Pr  {|GArrs  _r(k,,  v)  |  <  r\/(u.  v)  outside  main  beam} 

=  (l  -  e~Nr  J  .  (34) 

For  the  ARRS-P  and  ARRS-R,  the  variance  is  approximately 
independent  of  observation  angle,  so  we  must  observe  enough 
locations  that  the  joint  probability  will  be  representative.  For 
Lo’s  analysis  of  a  purely  random  array  the  number  of  locations 


observed  was  16  x  (array  area  in  square  wavelengths)  [4].  Ap¬ 
plying  this  to  the  ARRS-P,  we  obtain 

Pr  {|GARRS-p(tq  n)|  <  rV(w,.  v)  outside  main  beam} 

ss  (Erf(rv^)V))fa4a  .  (35) 

If  there  is  more  than  one  block  of  four  subarrays  in  an  ARRS-P, 
the  blocks  themselves  may  be  treated  as  subarrays  in  a  PARS. 
In  that  case,  the  probability  distribution  becomes 

Pr  {|GARRS-p(tb  n)|  <  rV(u,  v )  outside  main  beam} 

ss  ^Erf (rV2N)J  .  (36) 
Applying  the  same  method  to  the  ARRS-R 

Pr  { |Garrs-r(^;  v)  |  <  r,  V(n,  v)  outside  main  beam} 

/  PQa2 

=  (l  -  e~NPQr  )  .  (37) 

IV.  Array  Behavior 

The  goal  of  this  work  is  to  arrive  at  a  new  type  of  array  that 
duplicates  the  wideband  capabilities  of  purely  random  arrays 
while  improving  manufacturability  by  simplifying  the  geom¬ 
etry.  This  section  will  compare  how  well  the  three  new  types  of 
arrays  that  have  been  presented  meet  that  objective.  Section  Ill’s 
probabilistic  analyses  will  be  employed  to  compare  characteris¬ 
tics  such  as  the  number  of  elements  needed  to  meet  a  particular 
peak  SLL  requirement,  the  relative  peak  sidelobe  levels  achiev¬ 
able  by  arrays  with  the  same  number  of  elements  and  aperture 
size,  the  extent  to  which  division  into  subarrays  is  practical  for 
simplifying  the  geometry,  and  the  directivities  achievable  with 
each  type  of  array. 

A.  Probabilistic  Comparison  of  Array  Properties 

Figs.  5  and  6  show  the  peak  SLL  probability  distribution  for 
all  three  arrays  discussed  here  versus  r  for  arrays  of  four  and 
sixteen  subarrays,  respectively.  As  is  evident  in  these  figures, 
the  peak  sidelobe  levels  of  all  three  arrays  exhibit  threshold  be¬ 
havior  in  r.  There  is  a  threshold  where  the  probability  of  a  suc¬ 
cessful  array  increases  sharply  and  beyond  which  it  is  nearly  1 , 
increasing  very  slowly.  Similar  behavior  can  be  observed  in  the 
plots  of  peak  SLL  versus  N  (Figs.  7  and  8). 

This  is  the  foundation  of  the  wide  bandwidth  capability  of 
random  arrays;  if  N  and  r  are  chosen  to  be  above  threshold  over 
the  entire  frequency  range,  an  array  with  those  parameters  is  al¬ 
most  certain  to  have  acceptable  sidelobe  levels  throughout  that 
range.  Fig.  9  shows  the  variation  of  the  peak  SLL  probability 
distribution  with  subarray  dimension  a  (in  this  example  the  sub¬ 
arrays  are  all  square),  which  is  proportional  to  frequency.  In 
Fig.  9,  N  is  sufficient  for  the  ARRS-R  probability  to  remain  very 
close  to  1  over  the  entire  frequency  range  while  the  ARRS-P  and 
PARS  probabilities  fall  off  more  quickly. 

Comparing  the  peak  SLL  probability  distribution’s  depen¬ 
dence  on  r  and  N  provides  information  about  the  relative  side¬ 
lobe  levels  achievable  with  each  type  of  array.  For  an  array  of 
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Probability 


Fig.  5.  Peak  SLL  probability  distribution  versus  r  for  PARS  (solid  line), 
ARRS-P  (dashed  line),  and  ARRS-R  (dot-dashed  line).  Each  array  has  four 
20A  X  20A  subarrays  of  36  elements. 


Probability 


Fig.  6.  Peak  SLL  probability  distribution  versus  r  for  PARS  (solid  line), 
ARRS-P  (dashed  line),  and  ARRS-R  (dot-dashed  line).  Each  array  has  sixteen 
20A  X  20A  subarrays  of  36  elements. 


Probability 


Fig.  8.  Peak  SLL  probability  distribution  versus  N  for  PARS  (solid  line), 
ARRS-P  (dashed  line),  and  ARRS-R  (dot-dashed  line).  Each  array  has  sixteen 
20A  X  20A  subarrays,  and  r  is  0.15. 


Probability 


Fig.  9.  Peak  SLL  probability  distribution  versus  a  for  PARS  (solid  line), 
ARRS-P  (dashed  line),  and  ARRS-R  (dot-dashed  line).  Each  array  has  four 
subarrays,  N  =  750  and  r  =  0.1. 


Probability 


Fig.  7.  Peak  SLL  probability  distribution  versus  N  for  PARS  (solid  line), 
ARRS-P  (dashed  line),  and  ARRS-R  (dot-dashed  line).  Each  array  has  four 
20A  X  20A  subarrays,  and  r  is  0.15. 


four  subarrays,  the  probability  distributions  for  the  ARRS-P  and 
ARRS-R  are  almost  identical,  which  is  expected  because  they 
work  on  the  same  principle.  In  the  case  of  four  subarrays,  the 
ARRS-P  effectively  counts  as  a  single,  diagonally  symmetric 
random  array,  so  its  probability  distribution  is  that  of  a  sym¬ 
metric  array  with  47V  elements.  As  the  number  of  subarrays 
grows,  however,  it  will  lose  this  advantage.  The  probability  dis¬ 
tribution  for  the  PARS  has  a  considerably  higher  threshold  be¬ 
cause  the  variance  of  its  probability  distribution  depends  only 
on  the  number  of  elements  in  the  subarray.  Fig.  6  demonstrates 


what  happens  to  the  ARRS-P  when  there  are  more  than  four  sub¬ 
arrays.  Recall  (36),  the  probability  distribution  for  an  ARRS-P 
with  more  than  four  subarrays.  Like  the  PARS,  the  ARRS-P 
gains  no  further  advantage  in  SLL  as  the  number  of  subarrays 
increases.  The  ARRS-R,  on  the  other  hand,  gets  a  substantial 
improvement  in  the  threshold  of  its  probability  distribution  due 
to  the  increase  in  number  of  subarrays. 

These  figures  show  probability  distributions  for  all  three  array 
types  for  four-subarray  and  sixteen-subarray  cases.  Clearly,  as 
the  number  of  subarrays  increases,  the  sidelobe  levels  of  the 
ARRS-P  and  PARS  essentially  stop  improving,  while  that  of  the 
ARRS-R  can  achieve  lower  sidelobe  levels  with  fewer  elements 
per  subarray  as  more  subarrays  are  added. 

Figs.  10  and  1 1  plot  the  number  of  elements  required  per  sub¬ 
array  for  PARS,  ARRS-P,  and  ARRS-R  arrays  with  four  and 
sixteen  subarrays,  respectively.  Again,  the  abilities  of  the  PARS 
and  ARRS-P  to  achieve  low  sidelobe  levels  are  not  improved  by 
the  addition  of  more  subarrays.  The  ARRS-R,  on  the  other  hand, 
has  the  capability  to  reduce  the  number  of  elements  that  need  to 
be  randomly  distributed  while  maintaining  the  frequency  and 
probabilistic  characteristics  of  a  purely  random  array. 

B.  Directivity  Comparison 

Random  array  design  has  historically  been  focused  on 
achieving,  with  fewer  elements  than  traditional  designs,  the 
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Fig.  10.  N  required  for  85%  probability  of  r  =  0.1,  versus  a  for  an  array  of 
four  subarrays  (PARS:  solid  line,  ARRS-P:  dashed  line,  ARRS-R:  dot-dashed 
line). 


Directivity 


Fig.  12.  Directivity  (in  dB)  of  ARRS-P  (dashed)  and  ARRS-R  (dot-dashed) 
versus  their  total  number  of  elements,  compared  to  that  of  a  periodic  array  with 
40,000  elements,  held  constant.  All  three  arrays  are  square  and  have  constant 
aperture  area  =  10,  000A2;  the  periodic  array  has  A/ 2  spacing. 


N  Required 


Fig.  11.  N  required  for  85%  probability  of  r  =  0.1,  versus  a  for  an  array  of 
sixteen  subarrays  (PARS:  solid  line,  ARRS-P:  dashed  line,  ARRS-R:  dot-dashed 
line). 


Also,  if  the  main  beam  is  narrow,  the  normalized  power  radiated 
in  the  sidelobes  is 

Psl  «  27rr2SLLa„  (40) 

where  SLL,,,,  is  the  average  power  sidelobe  level,  which  is  equal 
to  r2vg.  Combining  these  three  equations,  directivity  is  approx¬ 
imated  as 


~  isin2(^)  +  iSLLa„' 

Since  the  ARRS-P  and  ARRS-R  are  the  most  promising  of 
the  arrays  presented,  we  will  compare  their  directivities  with 
that  of  a  uniformly  illuminated  rectangular  aperture  with  100% 
aperture  efficiency.  From  [15],  the  directivity  of  such  an  aper¬ 
ture  at  broadside  is 


narrow  beamwidths  and  low  sidelobe  levels  required  for  appli¬ 
cations  like  radio  astronomy.  As  a  result,  many  random  arrays 
in  the  past  have  had  low  element  density  and  consequently  low 
directivity.  However,  the  directivity  of  a  random  array  depends 
largely  on  the  element  density  its  designer  allows,  and  it  is 
possible  to  design  random  arrays  with  useful  directivities  for  a 
range  of  applications.  This  section  will  use  a  simple  directivity 
approximation  [14]  to  illustrate  the  effect  of  element  density 
on  directivity. 

The  directivity  of  an  antenna  can  be  defined  as 


D  = 


Air  r2 

I  Js  p{0,<f))r2sm(6)ddd<p 


(38) 


where  p(6,<f>)  is  the  normalized  radiated  power  [15].  The  ap¬ 
proximation  we  will  use  assumes  that,  for  an  antenna  with  a 
half-power  beamwidth  of  26\ ,  the  normalized  power  in  the  main 
beam  is 


Pmb 


4-7rr2 


sin2(0i). 


(39) 


-Aperiodic  —  4  7T  .4 .  (42) 

A  periodic  array  with  the  same  aperture  size  whose  element 
spacing  is  a  half  wavelength  has  the  same  directivity  as  this  aper¬ 
ture  [15]. 

Using  the  directivity  approximation  from  (41),  an  approxi¬ 
mate  beamwidth  [  15],  and  the  average  sidelobe  levels  from  (32) 
and  (33),  Fig.  12  compares  the  directivity  of  the  ARRS-P  and 
ARRS-R  with  that  of  the  uniformly  illuminated  aperture.  The 
ARRS-R  in  this  case  has  a  continuum  of  allowed  rotations.  All 
three  have  aperture  area  held  constant  at  10000  A2;  everything 
remains  constant  about  the  aperture  while  the  density  of  ele¬ 
ments  in  the  ARRS-P  and  ARRS-R  increases.  This  comparison 
is  not  entirely  realistic  since  it  is  assumed  that  the  density  can 
arbitrarily  increase;  in  reality  elements  have  a  finite  size  and  it 
becomes  more  likely  that  some  might  overlap  as  the  density  in¬ 
creases. 

Fig.  12  illustrates  that,  as  one  would  expect,  the  directivity 
of  an  ARRS-P  or  ARRS-R  goes  up  as  the  element  density  in¬ 
creases.  By  the  time  the  random  arrays  have  the  same  element 
density  as  the  half-wavelength  spaced  periodic  array,  their  di¬ 
rectivities  are  less  by  about  8.5  and  6  dB,  respectively.  When 
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the  element  density  is  much  lower,  for  instance  between  den¬ 
sities  of  one  element  per  five  square  wavelengths  and  one  ele¬ 
ment  per  one  square  wavelength,  the  directivity  is  between  21 
dB  (ARRS-P,  2000  total  elements)  and  1 1  dB  (ARRS-R,  10,000 
total  elements)  less  than  that  of  the  periodic  array  which  has 
four  to  twenty  times  more  elements.  This  trade-off  in  directivity 
should  be  considered  in  conjunction  with  the  greatly  increased 
bandwidth  and  reduced  number  of  elements  of  the  ARRS-P  and 
ARRS-R. 


V.  Design  Example 

It  may  be  illustrative  to  design  a  sample  array  of  each  type, 
given  some  arbitrary  specifications.  Suppose  that  an  array  is  de¬ 
sired  that  operates  from  2  to  10  GHz  and  has  a  peak  SLL  of  —20 
dB.  Because  the  analysis  of  these  arrays  depends  on  the  central 
limit  theorem,  the  number  of  elements  per  subarray  must  be  rea¬ 
sonably  large.  We  will  require  at  least  20.  Since  antennas  with 
wide  bandwidths  often  have  larger  dimensions,  we  will  restrict 
the  element  density  to  at  most  one  per  three  square  wavelengths 
to  prevent  overlapping.  A  nearest-neighbor  constraint  could  also 
be  imposed,  but  for  a  simple  example  the  density  requirement 
should  be  sufficient.  Because  smaller  subarrays  allow  greater 
simplification  of  feed  networks  and  the  overall  manufacturing 
process,  an  array  with  many  subarrays  of  a  small  number  of  el¬ 
ements  is  preferred  over  one  with  fewer  subarrays  and  more  el¬ 
ements  per  subarray. 

In  this  section,  plots  of  each  array  factor  will  be  presented 
for  only  the  highest  frequency  in  the  selected  operating  band. 
An  array  factor  can  be  considered  an  arbitrary  function  of  au 
and  av,  where  a,  the  subarray  dimension  in  wavelengths,  scales 
linearly  with  frequency.  For  any  such  function,  increasing  a  will 
scale  the  features  of  the  function  closer  together  in  u  and  v, 
but  the  features  themselves  will  be  otherwise  unchanged.  By 
this  reasoning,  a  plot  of  an  array  factor  at  the  highest  point  in 
the  operating  band  indicates  the  array’s  behavior  over  the  entire 
band. 

If  we  require  that  the  probability  of  success  be  85%,  to  re¬ 
duce  the  chances  of  our  having  to  test  more  than  one  array,  this 
set  of  requirements  results  in  a  PARS  with  1309  elements  per 
subarray.  Assuming  the  subarrays  are  square,  they  would  mea¬ 
sure  62.67  wavelengths  on  a  side  at  2  GHz,  or  about  9.40  m. 
Fig.  13  shows  the  radiation  pattern  of  this  array  at  10  GHz. 

In  the  case  of  an  ARRS-P,  the  elements  decrease  the  SLL 
most  efficiently  if  there  are  four  subarrays.  Predetermining  the 
number  of  subarrays  selects  whether  the  array  is  described  by 
(35)  or  by  (36).  Having  determined  that,  the  same  probability 
requirement  we  enforced  on  the  PARS  now  yields  an  ARRS-P 
with  664  elements  per  subarray.  These  subarrays  are  required 
to  be  square  because  of  the  rotation  pattern,  and  they  would 
measure  44.64  wavelengths  on  a  side  at  2  GHz,  or  about  6.70 
m.  The  entire  array  of  four  subarrays  would  also  be  square,  and 
measure  13.40  m  on  a  side.  The  total  number  of  elements  is 
2656.  Fig.  14  shows  the  radiation  pattern  of  this  array  at  10  GHz. 

Lastly,  the  ARRS-R  is  considered.  Since  the  ARRS-R  peak 
sidelobe  probability  distribution  can  be  improved  by  increasing 
the  total  number  of  elements,  rather  than  the  number  of  elements 
per  subarray,  we  would  like  these  subarrays  to  be  as  small  as 
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Fig.  13.  Array  factor  at  10  GHz  of  PARS  designed  using  parameters  from 
sample  design:  N  =  1309,  a  =  b  =  634.5Ai0  ghz>  4  subarrays,  normal 
element  distribution. 
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Fig.  14.  Array  factor  at  10  GHz  of  ARRS-P  designed  using  parameters  from 
the  sample  design:  N  =  664,  a  =  b  =  223.20Aio  ghz>  4  subarrays,  normal 
element  distribution. 


possible  because  that  furthers  our  goal  of  simplifying  the  geom¬ 
etry.  Assuming  there  are  20  elements  per  subarray,  the  resultant 
array  has  72  subarrays  and  a  total  of  1440  elements.  Because  we 
must  allow  arbitrary  rotation  amounts  without  overlapping,  the 
subarrays  are  assumed  to  be  circular  since  that  minimizes  the 
surface  area.  The  element  density  requirement  results  in  subar¬ 
rays  with  a  diameter  of  1.31  m.  This  distance  will  also  be  the 
spacing  a  of  the  subarray  grid.  If  the  array  is  a  rectangle  of  9 
subarrays  by  8  subarrays,  its  dimensions  are  10.5  by  11.8  m. 
Fig.  15  shows  the  radiation  pattern  of  this  array  at  10  GHz. 

Fig.  16  shows  the  main  beam  region  of  the  array  factor 
at  2  GHz  of  an  ARRS-R  with  properties  determined  by  the 
sample  design.  Superarray  grating  lobes  are  clearly  visible  near 
the  main  beam.  This  occurs  because  the  array  factor  of  the 
subarray  is  dominated  by  its  expected  value  in  that  region,  so 
the  subarray  patterns  are  still  too  similar  after  rotation  to  reduce 
the  grating  lobes  to  the  same  degree  that  they  are  reduced  in  the 
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Fig.  15.  Array  factor  at  10  GHz  of  ARRS-R  designed  using  parameters  from 
sample  design:  N  =  20 ,  a  =  43.70Aio  ghz>  72  subarrays,  normal  element 
distribution. 


array  factor  detail 
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Fig.  1 6.  Main  beam  region  of  array  factor  at  2  GHz  of  ARRS-R  designed  using 
parameters  from  sample  design:  N  =  20,  a  =  8.74A2ghz>  72  subarrays, 
normal  element  distribution. 


far  sidelobe  region.  It  is  important  to  note  that,  since  the  shape 
of  the  subarray  main  beam  depends  strongly  on  the  element 
distribution,  a  tapered  element  distribution  in  the  subarray  is 
likely  to  result  in  more  prominent  first  grating  lobes  because 
the  subarray  main  beam  will  be  wider.  Indeed,  if  one  compares 
Fig.  17,  in  which  the  subarrays  have  a  uniform  element  distri¬ 
bution,  with  Fig.  16,  where  the  element  distribution  is  normal, 
it  is  apparent  that  fewer  grating  lobes  are  visible  above  the 
random  sidelobe  variation. 

It  is  probably  not  possible  to  entirely  eliminate  these  grating 
lobes  in  this  type  of  array,  since  the  first  grating  lobe  loca¬ 
tion  for  the  superarray  pattern  is  very  near  the  predicted  half¬ 
power  beamwidth  of  a  subarray  with  uniform  element  distribu¬ 
tion  (Fig.  18).  However,  the  grating  lobes  have  significantly  less 
than  half  the  power  of  the  main  beam,  which  implies  that  the 
subarray  rotation  is  having  some  effect  on  the  SLL. 


array  factor  detail 


u 


Fig.  17.  Main  beam  region  of  array  factor  at  2  GHz  of  ARRS-R  designed  using 
parameters  from  sample  design:  N  =  20,  a  =  8.74A2ghz>  72  subarrays, 
uniform  element  distribution. 


angle 


Fig.  18.  Beamwidth  between  first  nulls  (dotted  line)  and  half-power  beamwidth 
(solid  line)  for  subarray  with  dimension  d,  compared  to  first  grating  lobe  loca¬ 
tion  of  superarray  with  grid  spacing  d  (dot-dashed  line).  Angles  in  radians. 


VI.  Conclusion 

This  investigation  has  yielded  two  useful  new  arrays:  the 
ARRS-P  arrangement  can  be  used  for  an  array  with  a  small 
number  of  subarrays,  and  the  ARRS-R  generalizes  the  rotation 
approach  to  an  arbitrary  number  of  subarrays.  The  arrays  were 
intended  to  update  the  random  array  concept  for  the  present-day 
demands  of  fabrication  while  preserving  its  bandwidth  capa¬ 
bility.  The  ARRS-R  has  accomplished  that  goal  quite  well, 
transforming  the  purely  random  array  into  many  small  identical 
subarrays  while  almost  duplicating  the  purely  random  array’s 
performance. 

Since  none  of  these  arrays  has  been  constructed,  one  por¬ 
tion  of  planned  future  work  is  an  experiment  much  like  the 
holey-plate  experiment  performed  by  Lo  [16].  To  reduce  the  size 
of  the  array  so  that  measurement  is  practically  feasible,  the  array 
will  be  scaled  down  and  the  experiment  performed  in  the  optical 
regime  [17],  [18].  Since  the  utility  of  the  ARRS-R  is  improved 
if  the  subarrays  are  allowed  to  have  fewer  elements,  one  future 
investigation  of  this  topic  could  focus  on  the  effect  of  small  N 
on  the  accuracy  of  this  analysis.  Also,  as  was  mentioned  in  the 
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introduction,  the  rotation  method  can  be  applied  to  other,  pos¬ 
sibly  optimized  types  of  aperiodic  arrays. 
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